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1 Introduction 

The exploration of the phase diagram of dense baryonic matter is an area of 
intense theoretical and experimental activity. Baryonic systems, from dilute 
neutron matter at low density to superconducting quark matter at high den- 
sity, exhibit an enormous variety of many-body effects. Despite its simplicity 
all these phenomena are ultimately described by the lagrangian of QCD. 

In practice it is usually very difficult to describe QCD many-body sys- 
tems directly in terms of the QCD lagrangian, and even in cases where this 
is possible it is often not the most convenient and most transparent descrip- 
tion. Instead, it is advantageous to employ an effective field theory (EFT) 
formulated in terms of the relevant degrees of freedom. EFTs also provide a 
unified description of physical systems involving very different length scales, 
such as Fermi liquids in nuclear and atomic physics, or non- Fermi liquid gauge 
theories involving colored quarks or charged electrons. 

In these lectures we shall discuss the many body physics of several effective 
field theories relevant to the structure of hadronic matter. We will concentrate 
on two regimes in the phase diagram. At low baryon density the relevant 
degrees of freedom are neutrons and protons, while at very high baryon density 
the degrees of freedom are quarks and gluons. These lectures do not provide 
an introduction to effective field theories (see 13 E]), nor an exhaustive 
treatment of many body physics (see 0JE1 Ej) or the physics of dense quark 
matter (see [HE])- 



2 Fermi liquids 

2.1 Effective field theory for non-relativistic fermions 

If the relevant momenta are small neutrons and protons can be described as 
point-like non-relativistic fermions interacting via local forces. Effective field 
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theories for nuclear systems have been studied extensively over the past couple 
of years |5ll91 llOUll| . If the typical momenta are on the order of the pion mass 
pions have to be included as explicit degrees of freedoms. For simplicity we 
will consider neutrons only and focus on momenta small compared to m n . 
The effective lagrangian is 

Co = tf + (V^) 2 + § [W) t (V'V 2 ^) + / l . C ]+..., (1) 

where to is the neutron mass, Cq and C2 arc dimcnsionful coupling constants, 

V=V — V is a Galilei invariant derivative, and . . . denotes interactions with 
more derivatives. We have only displayed terms that act in the s-channel. The 
coupling constant are determined by the neutron-neutron scattering ampli- 
tude. For non-relativistic scattering the amplitude is related to the scattering 
phase shift S by 

TO p COt — ip 

For small momenta the quantity p cot 6 can be expanded as a Taylor series in 
p. This expansion is called the effective range expansion 

1 1 00 

p CO t5 = - 1 -+ l -Y,r n p 2{n+1 \ (3) 

n=0 

where a is the scattering length, and ro is the effective range. The situation is 
simplest if the scattering length is small. In this case the scattering amplitude 
has a perturbative expansion in C{. At tree level 

c = ^, c 2 = c ^. (4) 

TO 2 

However, there are many systems of physical interest in which the scattering 
length is not small. This happens whenever there is a two-body bound state 
with a very small binding energy, or if the two-body system is very close to 
forming a bound state. For neutrons a nn = — 17 fm, much larger than typical 
strong interaction length scales. 

If the scattering length is large then loop diagrams with the leading order 
interaction Cq(i() ip) 2 bave to be resummed. The one-loop correction involves 
the loop integral 



L{E) 



d d+1 Q 1 



{2n) d + 1 (E/2 + q - q 2 /(2m) + ie)(E/2 - q - q 2 /{2m) + ie) 
d d q 1 



(2ir) d E — q 2 /m + ie 
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where E is the center-of-mass energy. We have regularized the integral by 
analytic continuation to d+1 dimensions. In order to define the theory we have 
to specify a subtraction scheme. Here, we will employ the modified minimal 
subtraction MS scheme. See f° r a discussion of different renormalization 
schemes. We get 

L(E) = —V-mE -ie = -— ip , (6) 

All 47T 

where p = V mE is the nucleon momentum in the center-of-momentum frame. 
It is now straightforward to sum all the bubble diagrams. The result is 

A = 9l ( 7 ) 

l + iCoimp/An)' { > 

Higher order corrections due to the C, terms (i > 2) can be treated perturba- 
tively. The bubble sum can now be matched to the effective range expansion. 
In the MS scheme the result is particularly simple since equ. © only contains 
the contribution from the unitarity cut. As a consequence, the result given in 
equ. J3J is not modified even if Co is summed to all orders. 

2.2 Dilute Fermi liquid 

The lagrangian given in equ. Q is invariant under the U(l) transformation 
ip — > e^ip. The U(l) symmetry implies that the fermion number 

N = I d 3 x^ijj (8) 



is conserved. As a consequence, it is meaningful to study a system of fcrmions 
at finite density p — N/V. We will do this in the grand-canonical formalism. 
We introduce a chemical potential /i conjugate to the fermion number N and 
study the partition function 

Z(^/3) = Tr\e-^ H -^] . (9) 

Here, H is the Hamiltonian associated with C and (3 = 1/T is the inverse 
temperature. The trace in equ. runs over all possible states of the system. 
The average number of particles for a given chemical potential /i and temper- 
ature T is given by (N) = T(d\ogZ) / (dp). At zero temperature the chemical 
potential is the energy required to add one particle to the system. 

There is a formal resemblance between the partition function equ. © and 
the quantum mechanical time evolution operator U = exp(— iHt). In order to 
write the partition function as a time evolution operator we have to identify 
/3 — > it and add the term —pN to the Hamiltonian. Using standard techniques 
we can write the time evolution operators as a path integral |13l I14| 

Z = J DipLhft exp J dr J d 3 xC E j ■ (10) 
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Fig. 1. Leading order Feynman diagrams for the ground state energy of a dilute 
gas of fermions interacting via a short range potential. 

Here, Ce is the euclidean lagrangian 

72 



(^-A*-^) ^ + ^(V- f V-) 2 + ... ■ (11) 



The fermion fields satisfy anti-periodic boundary conditions i/j((3) = —ip(0). 
Equation (jl 1|> is the starting point of the imaginary time formalism in thermal 
field theory. The chemical potential simply results in an extra term —fitp^ip 
in the lagrangian. From equ. we can easily read off the free fermion 
propagator 

S° aP (P) = . p2 ■ (I 2 ) 

IP4 + M — 2m 

where a, (3 are spin labels. We observe that the chemical potential simply 
shifts the four-component of the momentum. This implies that we have to 
carefully analyze the boundary conditions in the path integral in order to fix 
the pole prescription. The correct Minkowski space propagator is 

S%{P) = H , ; =6 a J g^Zgl + e <PF-P} \, (13) 

Po - e P + iosga(e p ) I po - e p + 10 p - t p - 10 J 

where e p = E p — /.i, E p = p 2 /(2m) and 6 — * + . The quantity pf — y/2mfj, is 
called the Fermi momentum. We will refer to the surface defined by the con- 
dition |p| = pf as the Fermi surface. The two terms in equ. (|13fl have a simple 
physical interpretation. At finite density and zero temperature all states with 
momenta below the Fermi momentum are occupied, while all states above the 
Fermi momentum are empty. The possible excitation of the system are parti- 
cles above the Fermi surface or holes below the Fermi surface, corresponding 
to the first and second term in equ. I|13|l . The particle density is given by 

p = <v^> = / ^-Aip) <? vo V 0+ = 2 / (03 (^ -p) = ( 14 ) 



Tadpole diagrams require an extra iS prescription which can be derived from a 
careful analysis of the path integral representation at fi 7^ 0. As a first simple 
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application we can compute the energy density as a function of the fermion 
density. For free fermions, we find 



We can also compute the corrections to the ground state energy due to the 
interaction (Co/2)(ijj^ijj) 2 . The first term is a two-loop diagram with one inser- 
tion of Co, see Fig.^ There are two possible contractions and the spin- factor 
of the diagram is {5 aa 5fjp — S a p5 a f3) = g{g — 1) where g = (2s + 1) is the 
degeneracy and s is the spin of the fermions. In the following we will always 
set g = 2. The diagram is proportional to the square of the density and we 
get 



We observe that the sum of the first two terms in the energy density can be 
written as 



which shows that the Co term is the first term in an expansion in ppa, suitable 
for a dilute, weakly interacting, Fermi gas. 

2.3 Higher order corrections 

The expansion in (pfQ>) was carried out to order (ppa) 2 by Huang, Lee and 
Yang Since then, the accuracy was pushed to 0((pi?a) 4 \og(ppa)), 

see for an EFT approach to this calculation. The 0{(pfo) 2 ) calculation 
involves a few new ingredients and we shall briefly outline the main steps. 
Consider the third diagram in Fig. 2] The contribution to the vacuum energy 
is 



-4 /!!/!!/ 1|^)% 2 )%3)% 1+92 - 5 3). as) 



We begin by performing two of the energy integrals using contour integration. 
The contours can be placed in such a way that the two poles correspond to 
two particles or two holes (but not a particle and a hole). This allows us to 
write 




(15) 




(16) 




(17) 




(19) 



where 9 q is the Pauli-blocking factor corresponding to a pair of holes 



e, 



9(p F - q x ) 6{p F - q 2 ) , 



(20) 
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and IIpp is the one-loop particle-particle scattering amplitude. Since qi : 2 are 
on-shell we can write n pp as a function of the center-of-mass and relative mo- 
menta P and k with q 12 = P/2±k. Note that because of Galilean invariance 
the vacuum scattering amplitude only depends on k. We find 

TT - f d3<1 ^ - ttvociu\ i m P F * U .\ \ 

Upp ~ J (2^)8 k 2 -q 2 + ie ~ pp [ ' + (2^ /pp(K ' S > ' 1 > 

where 6>+ = 9 (pp + qi) 9 (pp + q 2 ) is defined in analogy with equ. (|2*U|) . The 
first term on the RHS is the vacuum contribution and the second term is 
the medium contribution which depends on the scaled momenta k — k/pp 
and s = P/(2pp). The vacuum contribution is divergent and needs to be 
renormalized. In dimensional regularization 77pp C is purely imaginary and does 
not contribute to the vacuum energy. In other regularization schemes the 
vacuum contributions combines with the 0(Cq) graph to give the correct one- 
loop relation between Co and the scattering length. 

For s < 1 the in-medium scattering amplitude is given by 



fpp{n, s) = 1 + s + k log 



1 



1 + S + K 



1 — K 



2 



2s 



■log 



(1 + S) 2 - K 2 



1 



(22) 



The contribution to the energy density can now be determined by integrating 
equ. l|22[l over phase space. We find 



2 p F m I d 3 P d 3 k 



= ^3i (11 - 21 ° g(2)) ^ a)2 - (23) 
The fourth diagram in Fig. ^ involves a particle-hole pair with zero energy 
and the corresponding phase space factor vanishes 5 . 

The effective lagrangian can also be used to study many other properties 
of the system, such as corrections to the fermion propagator. Near the Fermi 
surface the propagator can be written as 

s = ZS aP / 24 ) 

Po -v F (\p\ -pp) +iSsgn(\p\ -pf) 1 

where Z is the wave function renormalization and Dp = pp /m* is the Fermi 
velocity. Z and m* can be worked out order by order in (pfci), see |H H%]. 
The leading order results are 

— = 1 - (1 " 71og(2)) {ppaf + ... (25) 

Z- 1 = l-^\og(2) (ppaf + ... (26) 

The main observation is that the structure of the propagator is unchanged 
even if interactions are taken into account. The low energy excitations are 
quasi-particles and holes, and near the Fermi surface the lifetime of a quasi- 
particle is infinite. This is the basis of Fermi liquid theory |19l I2()j . 
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Fig. 2. Leading order Feynman diagrams that contribute to the photon polarization 
function in a non-relativistic Fermi liquid. The tadpole diagram shown in the right 
panel only appears in the spatial part of the polarization tensor. 



2.4 Screening and damping 

An important aspect of a dilute Fermi gas of charged particles is the response 
to an external electromagnetic field. We consider a system in which the total 
charge is neutralized by a homogeneous background (such as positive ions in 
a metal). The response to an electric field is governed by the gauge coupling 
eAQip^ip. The medium correction to the photon propagator is determined by 
the polarization function 

n o(q) =e 2 [ d^xe-^^mo^^x)). (27) 



The one-loop contribution is given by 



#oo (q) = -ie 



d A p 1 



(27r) 4 q + p - e p+q + i<Ssgn(e p+g ) p - e p + i5sga(e p ) ' 

(28) 

Performing the po integral using contour integration we find 

(29) 

where we have introduced the Fermi distribution function n p = 0{pf — p)- 
We observe that in the limit q — > the polarization function only receives 
contributions from particle-hole pairs that are very close to the Fermi surface. 
On the other hand, the energy denominator diverges in this limit because the 
photon can excite particle-hole pairs with arbitrarily small energy. These two 
effects combine to give a finite contribution 

W-0.,~Q-*f ,30, 

which is proportional to the density of states on the Fermi surface. Equ. lj3U|) 
implies that the static photon propagator in the limit q — > is modified 
according to l/q 2 — > 1/(<Z 2 +m 2 D ), where 



e 2 
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is called the Debye mass. The factor N = (pFm)/(2TT 2 ) is equal to the density 
of states on the Fermi surface. In a relativistic theory we find the same result 
as in equ. (|31[1 with the density of states replaced by the correct relativistic 
expression N = (ppEp)/ (2tt 2 ). The Coulomb potential is modified as 

p-r/rn 

V(r) = -e- , (32) 

r 

where m = l/mo is called the Debye screening length. The physics of screen- 
ing is very easy to understand. A test charge can polarize virtual particle-hole 
pairs that act to shield the charge. 

In the same fashion we can study the response to an external vector po- 
tential A. The coupling of a non-relativistic fermion to the vector potential is 
determined in the usual way by replacing p — ► p + eA. Since the kinetic en- 
ergy operator is quadratic in the momentum we find a linear and a quadratic 
coupling of the vector potential. The one-loop diagrams that contribute to the 
polarization tensor are shown in Fig. [2] In the limit of small external momenta 
we find 

n ij(v) = -e 2 m 2 D / — I — — h -v 2 5 lJ \ , (33) 

J ™ ( q -v ■ q 3 J 

where v = p/m is the Fermi velocity. In the limit q — the polarization 
tensor vanishes. There is no screening of static magnetic fields. For non-zero 
go the trace of the polarization tensor is given by 

n u (q) = m D —\og —— — . 34) 
2g \go + vqJ 

The result has an imaginary part for vq > q$. This phenomenon is known 
as Landau damping. The physical mechanism is that the photon is loosing 
energy as it scatters of the electrons in the Fermi liquid, see j2J f° r a detailed 
discussion in the context of kinetic theory. 



3 Superconductivity 
3.1 BCS instability 

One of the most remarkable phenomena that take place in many body systems 
is superconductivity. Superconductivity is related to an instability of the Fermi 
surface in the presence of attractive interactions between fermions. Let us 
consider fermion- fermion scattering in the simple model introduced in Sect. [21 
At leading order the scattering amplitude is given by 

r a f3ys(pi,P2,P3,P4) = Co (5 al S/3S ~ SasSfSy) . (35) 

At next-to-leading order we find the corrections shown in Fig. [3] A detailed 
discussion of the role of these corrections can be found in 0J|22in]- The BCS 
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Fig. 3. Second order diagrams that contribute to particle-particle scattering. The 
three diagrams are known as the ZS (zero sound), ZS' and BCS (Bardeen-Cooper- 
Schrieffer) contribution. 



diagram is special, because in the case of a spherical Fermi surface it can lead 
to an instability in weak coupling. The main point is that if the incoming 
momenta satisfy p\ ~ —P2 then there are no kinematic restrictions on the 
loop momenta. As a consequence, all back-to-back pairs can mix and there is 
an instability even in weak coupling. 

For pi = —p2 and E\ = E2 = E the BCS diagram is given by 



1 



d 4 q 

(2tt) 4 E + q - e q + iSsga(e q ) 



E - q a - e q + i8sgn(e q ) 



(36) 



As E — > the loop integral develops an infrared divergence. This divergence 
comes from momenta near the Fermi surface and we can approximate d 3 q ~ 
p 2 F dl with I = \q\ — pp. The scattering amplitude is proportional to 



Eo 
E 



(37) 



where Eq is an ultraviolet cutoff. The logarithmic divergence can also be seen 
by expanding equ. (|22|l around s — and n = 1. The term in the curly 
brackets can be interpreted as an effective energy dependent coupling. The 
coupling constant satisfies the renormalization group equation 122| 



dE 



with the solution 



C (E) 



Co(E ) 



l + NC {E )log(E /Ey 



(38) 



(39) 



where N — (p^m)/(27r 2 ) is the density of states. Equ. J3^|l shows that there 
are two possible scenarios. If the initial coupling is repulsive, Cq(Eq) > 0, 
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Fig. 4. Gap equation for the superfluid gap in a theory with short range interac- 
tions. 

then the renormalization group evolution will drive the effective coupling to 
zero and the Fermi liquid is stable. If, on the other hand, the initial coupling 
is attractive, Cq(Eq) < 0, then the effective coupling grows and reaches a 
Landau pole at 



At the Landau pole the Fermi liquid description has to break down. The renor- 
malization group equation does not determine what happens at this point, 
but it seems natural to assume that the strong attractive interaction will 
lead to the formation of a fermion pair condensate. The fermion condensate 
{e al3 ipa' t p0) signals the breakdown of the U(l) symmetry and leads to a gap 
A in the single particle spectrum. 

The scale of the gap is determined by the position of the Landau pole, 
A ~ E cr it . A more quantitative estimate of the gap can be obtained in the 
mean field approximation. In the path integral formulation the mean field ap- 
proximation is most easily introduced using the Hubbard-Stratonovich trick. 
For this purpose we first rewrite the four-fermion interaction as 



where we have used the Fierz identity 25 a f 3 8~<P = 5 a P6' fl3 + (tr) a f>(tr)iP . Note 
that the second term in equ. (j41(l vanishes because (<72 <t) is a symmetric 
matrix. We now introduce a factor of unity into the path integral 




(40) 



-fi^Hf = -f {OW f )(^2V0 + (vW^XVxr^V)} > (4i) 




(42) 



where we assume that Cq < 0. We can eliminate the four-fermion term in the 
lagrangian by a shift in the integration variable A. The action is now quadratic 
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in the fermion fields, but it involves a Majorana mass term i/jo^A^-I- h.c. The 
Majorana mass terms can be handled using the Nambu-Gorkov method. We 
introduce the bispinor W — (ip,^ ^2) and write the fermionic action as 

5 = 1 A V- (43) 

Since the fermion action is quadratic we can integrate the fermion out and 
obtain the effective lagrangian 

£=iTr[log(G - 1 G)]+-L|Z\| 2 , (44) 

Z Oo 

where G is the fermion propagator 

G(p)= a I lA J P °T P A * )■ ( 45 ) 
P - £ p- \ A \ V A Po-tpJ 

The diagonal and off-diagonal components of G(p) are sometimes referred to 
as normal and anomalous propagators. Note that we have not yet made any 
approximation. We have converted the fermionic path integral to a bosonic 
one, albeit with a very non-local action. The mean field approximation corre- 
sponds to evaluating the bosonic path integral using the saddle point method. 
Physically, this approximation means that the order parameter does not fluc- 
tuate. Formally, the mean field approximation can be justified in the large N 
limit, where N is the number of fermion fields. The saddle point equation for 
A gives the gap equation 

Performing the po integration we find 

1= mfijL 1 (47) 



Since e p = E p — fi the integral in equ. (|47|l has an infrared divergence on the 
Fermi surface \p\ ~ pf- As a result, the gap equation has a non-trivial solution 
even if the coupling is arbitrarily small. We can estimate the size of the gap 
as we did earlier by writing d 3 p ~ ppdl and introducing a cutoff A for the 
integral over /. We find A = 2Aexp(—l/(N\Co\))- In order to obtain a more 
accurate result we compute the RHS of equ. I|47|l without the approximation 
d 3 q ~ ppdl. We use dimensional regularization and 

dz - Z ° * ( 1+3 ?)"l 2 Pct ( 1 V (48) 

y/(z- l) 2 +x 2 sm(7ra) v > \ ^IT^ 2 J 

The dimensionally regularized gap equation is |231 12*1| 
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l=T^4(l + ^) a/9 Paf~™), (49) 

sm(Tra) \ \/l + x 2 J 

where 2A = Comp F ~ 2 {2^/ (2n) d is a dimensionless coupling constant, fi^ is the 
surface area of the <i-dimensional unit ball and x = AjEp is the dimensionless 
gap. P a (z) is the Legendre function of order a and a — (d—2)/2. Dimensional 
regularization sets the power divergence in equ. (|4*T|l to zero. As a result, we 
can set d = 3 and Co = 47ra/m in equ. l)15). If the gap is small, iCl, equ. if^l 
can be solved using the asymptotic behavior of the Legendre function P a {z) 
near the logarithmic singularity at z = — 1, 

P Q ( z ) ~ sm ^ ) hog f i±£ J + 2 7 + 2^(a + 1) + TT cot(aTr)) . (50) 

We find 

^"^"pf-rTi)- (51) 

e 2 V 2 PF|a|y 

The term in the exponent represents the leading term in an expansion in 
Pi?|a|, see Fig. This means that in order to determine the pre-exponent 
in equ. I|51|) we have to solve the gap equation at next-to-leading order. The 
contribution from the second diagram in Fig.^t was first computed by Gorkov 
and Melik-Barkhudarov . The second order graph screens the leading order 
particle-particle scattering amplitude and suppresses the s-wave gap by a 
factor (4c) 1 / 3 - 2.2 

For neutron matter the scattering length is large, a = —18.8 fm, and 
equ. 151(1 is not very useful, except at very small density. At moderate den- 
sity a rough estimate of the gap can be obtained by replacing l/(ppa) with 
cot(<5(fcp)), where S(k) is the s-wave phase shift. This estimate gives neutron 
gaps on the order of 1 MeV at nuclear matter density. 



3.2 Superfluidity 

Pairing leads to important physical effects. If the fermions are charged, pairing 
causes superconductivity. If the fermions are neutral, pairing leads to super- 
fluidity. We first discuss superfluidity. The superfluid order parameter (i^ip) 
breaks the U(l) symmetry and leads to the appearance of a Goldstone boson. 
The Goldstone boson field is defined as the phase of the order parameter 

(iH)) = \{iH>)\4*e. (52) 

In the following we shall construct an effective lagrangian for the Goldstone 
field if. The U{1) symmetry implies that the lagrangian can only depend on 
derivatives of if. The simplest possibility is 



£ = f((d f) 2 ~v 2 (d l f) 2 + ...), 



(53) 
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where vf is the coupling of d^ip to the U(l) current and v is the Goldstone bo- 
son velocity. This Lagrangian correctly describes the propagation of Goldstone 
modes and the coupling to external currents, but it does not respect Galilean 
invariance, and it does not describe the interaction between Goldstone modes 
|26l 1271 . Under Galilean transformations the fermion field transforms as 

tp(t, x) -> xp'(t, x) = e lmv - x xP(t, x - vt). (54) 

This implies that tp transforms as ip(t, x) — > <p(t,x — vt) + mv ■ x. We also 
observe that the chemical potential enters the microscopic theory like the time 
component of a U(l) gauge field. We can impose the constraints of Galilei 
invariance and £7(1) symmetry by constructing an effective lagrangian that 
only depends on the variable 

X = „-d tp-¥f£. (55) 

In the following it will be useful to consider a low energy expansion in which 
dof, d{ip are O(l) but higher derivatives didjtp, etc. are suppressed. In this 
case the leading order lagrangian contains arbitrary powers of X, but terms 
with derivatives of X are suppressed. The functional form of C(X) can be 
determined using the following simple argument. For constant fields ip = const 
the lagrangian C(X) — C(u) is equal to minus the thermodynamic potential 
Q. Since Q = —P, where P is the pressure, we conclude that C(X) — P(X). 

As an example consider superfluidity in a weakly coupled Fermi gas. 
At leading order the equation of state is that of a free Fermi gas, P = 
m 3 / 2 (2/i) 5 / 2 /(157r 2 ). The effective lagrangian is given by 

C = ^^{^ d0V -—) ■ (56) 

We can determine the Goldstone boson propagator as well as Goldstone boson 
interactions by expanding this result in powers of dop and diip. There are 
some predictions that are independent of the equation of state. Consider the 
effective theory at second order in (dtp), 

C = P(p) - ndotp + \^(do^) 2 - ^(^p) 2 + ..., (57) 

where we have used n = (dP)/(dfj). The Goldstone boson velocity is given by 

« 2 = ^ = ^. (58) 
m on op 

where p = nm denotes the mass density. We observe that the Goldstone boson 
velocity is given by the same formula as the speed of sound in a normal fluid. 
In a weakly interacting Fermi gas v 2 = Vp/3. 
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It is also instructive to study the relation to fluid dynamics in more detail. 
The equation of motion for the field (p is given by 

d n+ — V(nV<p) = 0, (59) 
to 

where we have defined n — P'(X). Equ. (|59|l is the continuity equation for 
the current j M = n(l,v s ) where we have identified the fluid velocity 

v s = ^. (60) 

TO 

In the hydrodynamic description the independent variables are n and v s . We 
can derive a second equation by using the identity dP = ndfi. We get 

d v s + ^-Vv 2 s = -i-V/j. (61) 

2 TO 

This is the Euler equation for non- viscous, irrotational fluid. The fact that the 
flow is irrotational follows from the definition of the velocity as the gradient 
of tp. We conclude that the low energy effective lagrangian is equivalent to 
superfluid hydrodynamics. 



3.3 Landau- Ginzburg theory 

In this section we shall study the properties of a superconductor in more 
detail. Superconductors are characterized by the fact that the U(l) symme- 
try is gauged. The order parameter <P — (€ a ^tp a tp8) breaks U(l) invariance. 
Consider a gauge transformation 

A^A^ + d^A. (62) 

The order parameter transforms as 

<P -> cxp(2ze/l)<Z>. (63) 

The breaking of gauge invariance is responsible for most of the unusual prop- 
erties of superconductors |2H1 EO] . This can be seen by constructing the low 
energy effective action of a superconductor. For this purpose we write the 
order parameter in terms of its modulus and phase 

&(x) = exp(2ie(f)(x))<P(x). (64) 



The field 4> corresponds to the Goldstone mode. Under a gauge transformation 
4>(x) — » 4>(x) + A(x). Gauge invariance restricts the form of the effective 
Lagrange function as 



(65) 
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There is a large amount of information we can extract even without knowing 
the explicit form of L s . Stability implies that = d^cf) corresponds to a 
minimum of the energy. This means that up to boundary effects the gauge 
potential is a total divergence and that the magnetic field has to vanish. This 
phenomenon is known as the Meissner effect. 

Equ. H65J1 also implies that a superconductor has zero resistance. The equa- 
tions of motion relate the time dependence of the Goldstone boson field to 
the potential, 

4>{x) = -V(x). (66) 

The electric current is related to the gradient of the Goldstone boson field. 
Equ. H66J) shows that the time dependence of the current is proportional to 
the gradient of the potential. In order to obtain a static current the gradient 
of the potential has to vanish throughout the sample, and the resistance is 
zero. 

In order to study the properties of a superconductor in more detail we have 
to specify L s . For this purpose we assume that the system is time- independent, 
that the spatial gradients are small, and that the order parameter is small. In 
this case we can write 

L s = J d 3 x |-i|(V-2 i eA)<?| 2 + im2 f (^^) 2 -i g (^<P) 4 + ...|, 

(67) 

where ran and g are unknown parameters that depend on the temperature. 
Equ. H67fl is known as the Landau-Ginzburg effective action. Strictly speaking, 
the assumption that the order parameter is small can only be justified in the 
vicinity of a second order phase transition. Nevertheless, the Landau-Ginzburg 
description is instructive even in the regime where t — (T—T c ) /T c is not small. 
It is useful to decompose <P — p exp(2ie<f)) . For constant fields the effective 
potential, 

V{p) = -\rnlp 2 + \gp\ (68) 

is independent of (j). The minimum is at p\ = m H / g and the energy den- 
sity at the minimum is given by £ — —ra H /{Ag). This shows that the two 
parameters ttlh and g can be related to the expectation value of <S> and the 
condensation energy. We also observe that the phase transition is character- 
ized by m#(T c ) = 0. 

In terms of 4> & n d p the Landau-Ginzburg action is given by 

L s = J d 3 x j- 2e y (V0 - Af + l -m H p 2 - \gp* - \ (V P ) 2 | . (69) 

The equations of motion for A and p are given by 

VxB = 4e 2 p 2 (V0 - A) , (70) 
V 2 p = —m 2 H p 2 + gp 3 + 4e 2 p (V0 — A) 2 . (71) 
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Fig. 5. Leading order Feynman diagrams that contribute to the photon polarization 
function in a superconducting Fermi gas. The figure does not show the tadpole 
diagram. 

Equ. I|7UI) implies that V 2 £? = ~4e 2 p 2 B. This means that an external mag- 
netic field B decays over a characteristic distance A = l/(2ep). Equ. H71|) 
gives V 2 /? = —m 2 H p+. ... As a consequence, variations in the order parameter 
relax over a length scale given by £ = 1/to#. The two parameters A and £ are 
known as the penetration depth and the coherence length. 

The relative size of A and £ has important consequences for the properties 
of superconductors. In a type II superconductor £ < A. In this case magnetic 
flux can penetrate the system in the form of vortex lines. At the core of a 
vortex the order parameter vanishes, p = 0. In a type II material the core is 
much smaller than the region over which the magnetic field goes to zero. The 
magnetic flux is given by 



and quantized in units of nh/e. In a type II superconductor magnetic vortices 
repel each other and form a regular lattice known as the Abrikosov lattice. 
In a type I material, on the other hand, vortices are not stable and magnetic 
fields can only penetrate the sample if superconductivity is destroyed. 

3.4 Microscopic calculation of the screening mass 

In this section we shall study screening of gauge fields in a superconductor 
from a more microscopic point of view. The calculation is analogous to the one 
discussed in Sect. 12.41 The difference is that the propagators contain the gap, 
and that there is an extra trace over Nambu-Gorkov indices, see Fig. [5] The 
polarization functions contains normal contributions proportional to GnGn 
and G22G22 as well as anomalous terms proportional to G12G21, where Gij is 
the Nambu-Gorkov propagator give in equ. (|45J) . The sum of the normal and 
anomalous diagrams is given by 



The integral over po can be done by contour integration. The two terms in 
equ. {EH give equal contributions. We find 




(72) 




(73) 
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g °°( g = 0) = e V(^ ( £ 3 + ^)3/ 2 - (^) 

This integral is dominated by very small energies |e p | = \E P — fx\ ~ and we 
can approximate e p = vp{p — pf)- We find 

7Too(g = 0) = e 2 ^J, (75) 

which is identical to the result in the normal phase. There are a number of 
subtleties that are worth commenting on. First we note that the polarization 
function in the superfluid phase is analytic in the external momenta and 
we can set qo = q = from the beginning. We also note that the normal 
contribution is formally ultraviolet divergent. The correct prescription to deal 
with this divergence is to perform the po integral first '4: . Finally we observe 
that while the screening masses in the normal and superfluid phase are the 
same, only half of the result in the superfluid phase is contributed by the 
normal term. 

The calculation of the electric polarization function is easily generalized 
to the magnetic case. There are three diagrams. The first is the tadpole con- 
tribution discussed in Sect. 12.41 This contribution is proportional to the total 
density and is the same in the normal and superfluid phase. The normal and 
anomalous one-loop diagrams are similar to the electric case, but the cou- 
pling e 2 is replaced by e 2 ViVj in the normal contribution and e 2 Vi(—Vj) in 
the anomalous term. As a result the two terms cancel and the polarization 
function is given by the tadpole term 

77 y -(^0H- e 2 4%^?. (76) 

We find that there is a non-zero magnetic screening mass in the superfluid 
phase, and that the Meissner mass is controlled not by the gap, but by the den- 
sity of states on the Fermi surface. This does not contradict the fact that the 
magnetic screening mass goes to zero as A — ► 0. We find that the photon mass 
term has the structure m^A 2 , — vpA 2 /?>). This result can also be obtained by 
gauging the effective Lagrangian for the Goldstone boson, equ. together 
with the result v 2 = Vp /3 for the speed of sound in a weakly interacting Fermi 
gas. 



4 Strongly interacting fermions 

Up to this point we have concentrated on weakly coupled many body systems. 
In this section we shall consider a cold, dilute gas of fermionic atoms in which 
the scattering length a of the atoms can be changed continuously. This system 
can be realized experimentally using Feshbach resonances, see for a review. 
A small negative scattering length corresponds to a weak attractive interaction 
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Fig. 6. The left panel shows the scattering length of 40 K Atoms as a function of 
the magnetic field near a Feshbach resonance, from Regal (2005). The right panel 
shows the nucleon-nucleon scattering length in the So channel as a function of the 
pion mass. The scatter plot indicates the uncertainty due to higher order terms in 
the chiral effective lagrangian. Figure from Beane & Savage (2003). 



between the atoms. This case is known as the BCS limit. As the strength of the 
interaction increases the scattering length becomes larger. It diverges at the 
point where a bound state is formed. The point a = oo is called the unitarity 
limit, since the scattering cross section saturates the s-wave unitarity bound 
a = in/k 2 . On the other side of the resonance the scattering length is positive. 
In the BEC limit the interaction is strongly attractive and the fermions form 
deeply bound molecules. 

A dilute gas of fermions in the unitarity limit is a strongly coupled quantum 
liquid that exhibits many interesting properties. One interesting feature is 
universality. We are interested in the limit (kpa) — > oo and (kpr) — > 0, where 
Uf is the Fermi momentum, a is the scattering length and r is the effective 
range. From dimensional analysis it is clear that the energy per particle at 
zero temperature has to be proportional to energy per particle of a free Fermi 
gas at the same density 




The constant £ is universal, i. e. independent of the details of the system. 
Similar universal constants govern the magnitude of the gap in units of the 
Fermi energy and the equation of state at finite temperature. 

Universal behavior in the unitarity limit is relevant to the physics of dilute 
neutron matter. The neutron-neutron scattering length is a nn — —18 fm and 
the effective range is r nn = 2.8 fm. This means that there is a range of densities 
for which the inter-particle spacing is large compared to the effective range 
but small compared to the scattering length. It is interesting to note that the 
neutron scattering length depends on the quark masses in a way that is very 
similar to the dependence of atomic scattering lengths on the magnetic field 
near a Feshbach resonance [22] > see Fig. 
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Fig. 7. Lattice results for the energy per particle of a dilute Fermi gas from Lee & 
Schafer (2005). We show the energy per particle in units of 3Ef/5 as a function of 
temperature in units of Tf- 



4.1 Numerical Calculations 

The calculation of the dimensionless quantity £ is a non-perturbative problem. 
In this section we shall tackle this problem using a combination of effective 
field theory and lattice field theory methods. We will study an analytical 
approach in the next section. We first observe that in the low density limit 
the details of the interaction are not important. The physics of the unitarity 
limit is captured by an effective lagrangian of point-like fermions interacting 
via a short-range interaction. The lagrangian is 

c = ^ v^ + SV'T^ 2 ' (78) 

as in Equ. . The usual strategy for dealing with the four-fermion interaction 
is to use a Hubbard-Stratonovich transformation as in Sect. 13. ll The partition 
function can be written as |33j 

Z = J DsDcDc exp [-S] , (79) 

where s is the Hubbard-Stratonovich field and c is a Grassmann field. S is a 
discretized euclidean action 

S = [e' fiat c*(n)ci(n + 6) - e v - c » a ' s W+¥(i _ &h)c*{n)a{n) 

n.i 

-hJ2 [tf + Is) + c*{n)c t (n - f s )] + ~ ^ s 2 (n). (80) 

n,l s ,i n 
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Here i labels spin and n labels lattice sites. Spatial and temporal unit vectors 
are denoted by l s and 6, respectively. The temporal and spatial lattice spacings 
are b T and b. The dimensionless chemical potential is given by fi = \xb T . 
We define at as the ratio of the temporal and spatial lattice spacings and 
h = at /(2m). Note that for Cq < the action is real and standard Monte 
Carlo simulations are possible. 

The four-fermion coupling is fixed by computing the sum of all particle- 
particle bubbles as in Sect. 12. ll but with the elementary loop function regu- 
larized on the lattice. Schematically, 

mil v ^ 1 

I^ = q; + 22-v (81) 
p 

where the sum runs over discrete momenta on the lattice and E p is the lattice 
dispersion relation. A detailed discussion of the lattice regularized scattering 
amplitude can be found in |341 1351 I3"3"| . For a given scattering length a the 
four-fermion coupling is a function of the lattice spacing. The continuum limit 
correspond to taking the temporal and spatial lattice spacings b T , b to zero 

b T (i -> 0, bn 1/3 -> 0, (82) 

keeping an 1 / 3 fixed. Here, [i is the chemical potential and n is the density. 
Numerical results in the unitarity limit are shown in Fig. [7] From these simula- 
tions we concluded that £ = (0.09—0.42). Lee performed canonical simulations 
at T = and obtained £ = 0.25. Green Function Monte Carlo calcula- 
tions give [37| £ = 0.44, and finite temperature lattice simulations have been 
extrapolated to T = to yield similar results [381 EH| ■ 

4.2 Epsilon Expansion 

It is also desirable to find a systematic analytical approach to the dilute Fermi 
liquid in the unitarity limit. Various possibilities have been considered, such 
as an expansion in the number of fermion species |401 141) or the number of 
spatial dimensions 02 SB| ■ Nussinov & Nussinov observed that the fermion 
many body system in the unitarity limit reduces to a free Fermi gas near d = 2 
spatial dimensions, and to a free Bose gas near d = 4 |44j . Their argument was 
based on the behavior of the two-body wave function as the binding energy 
goes to zero. For d = 2 it is well known that the limit of zero binding energy 
corresponds to an arbitrarily weak potential. In d = 4 the two-body wave 
function at a = oo has a 1/r 2 behavior and the normalization is concentrated 
near the origin. This suggests the many body system is equivalent to a gas of 
non- interacting bosons. 

A systematic expansion based on the observation of Nussinov & Nussinov 
was studied by Nishida and Son |451 . In this section we shall explain their 
approach. We begin by restating the argument of Nussinov & Nussinov in the 
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effective field theory language. In dimensional regularization a — > oo corre- 
sponds to Co — > oo. The fermion-fermion scattering amplitude (see equ. EJ is 
given by 



•A(p ,p) 



-in 



J/2 



r 1 



(83) 



where (5 — * 0+. As a function of d the Gamma function has poles at d = 
2,4,... and the scattering amplitude vanishes at these points. Near d = 2 the 
scattering amplitude is energy and momentum independent. For d = 4 — e we 
find 



•A(p ,p) 



8ir A 



0(e 2 ) . 



(84) 



to 2 po — Ep/2 + iS 

We observe that at leading order in e the scattering amplitude looks like 
the propagator of a boson with mass 2m. The boson-fermion coupling is 
g 2 = (8ir 2 e)/m 2 and vanishes as e — ► 0. This suggests that we can set up 
a perturbative expansion involving fermions of mass m weakly coupled to 
bosons of mass 2m. In the unitarity limit the Hubbard-Stratonovich trans- 
formed lagrangian reads 



id + a 3 



72 1 



2m 



& + fi&^a 3 ]P + {W^a + W4> + h. 



(85) 



where 'P = (^|,"0|) T is a two-component Nambu-Gorkov field, <Xj are Pauli 
matrices acting in the Nambu-Gorkov space and a± = {a\ ± ia-i)j2. In the 
superfluid phase <f> acquires an expectation value. We write 



fJ 



\/87r 2 e / m0o V^ 4 



2tt 



(86) 



where 0o = (<fi)- The scale M 2 = m^o/(27r) was introduced in order to have 
a correctly normalized boson field. The scale parameter is arbitrary, but this 
particular choice simplifies some of the loop integrals. In order to get a well 
defined perturbative expansion we add and subtract a kinetic term for the 
boson field to the lagrangian. We include the kinetic term in the free part of 
the lagrangian 



2m 



00 (o"H 



!? + <p* (id + ^) ip . (87) 



The interacting part is 



Ci = g (^a + <pLp + h.c) + ^ - ip^ ( id 



4m 



(88) 



Note that the interacting part generates self energy corrections to the boson 
propagator which, by virtue of equ. H84JI. cancel against the kinetic term of 
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Fig. 8. Leading order contributions to the effective potential in the e expansion. 
Solid lines are fermion propagators, dashed lines are boson propagators, and the 
cross is an insertion of the chemical potential. 



boson field. We have also included the chemical potential term in £/. This is 
motivated by the fact that near d = 4 the system reduces to a non-interacting 
Bose gas and /i — ► 0. We will count fj, as a quantity of O(e). 

The Feynman rules are quite simple. The fermion and boson propagators 

are 



G(po,p) = 
D(p ,p) = 



Po ~ ~ 

i 

Po - E p /2 ' 



Po + E p -0o 
-0o Pa - E p 



(89) 
(90) 



and the fermion-boson vertices are iga . Insertions of the chemical potential 
are i/jta^. Both g 2 and fi are corrections of order e. In order to verify that 
the e expansion is well defined we have to check that higher order diagrams 
do not generate powers of 1/e. Studying the superficial degree of divergence 
of diagrams involving the propagators given in equ. I|89() one can show that 
there is only a finite number of one-loop diagrams that generate 1/e terms. 

The leading order diagrams that contribute to the effective potential are 
shown in Fig. EI The first diagram is the free fermion loop which is 0(1). The 
second diagram is the fi insertion which is 0(1) because the loop diagram is 
divergent in d = 4. The two-loop diagram is O(e) because of the factor of g 2 
from the vertices. The free fermion loop diagram is 



V 



(91) 



The integral can be computed analytically. Expanding to first order in e = 4— d 
we get 



Vo = ^ 



1 



7 - 3( 7 + log(2)) 
6 



2tt 



d/2 



(92) 



The /i insertion is given by 
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Vi=ft 



d d p 



En 



Again, the integral can be computed analytically. The result is 



(93) 



l-2( 7 -Iog(2)) 



™</>0 

2tt 



d/2 



(94) 



Nishida and Son also computed the two-loop contribution shown in Fig. |SJ 
The result is 

' X \ d / 2 

m<po\ 



V 2 = -Ce 



2tt 



(95) 



where C ~ 0.14424. We can now determine the minimum of the effective 
potential. We find 



o = ^ [l + (3C-l + log(2)) e + 0(e 2 )] 



(96) 



The value of V — Vq + V\ + V2 at </>o determines the pressure and n = dPjd[i 
gives the density. We find 



l-i(2 7 -l-log(2)) + 0(e 2 ) 



mcf)o 
2tt 



d/2 



For comparison, the density of a free Fermi gas in d dimensions is 

_ 2 kp 



(47r) d / 2 r (1 + |) 



(97) 



(98) 



This equation determines the relation between ep = k F / (2m) and the density. 
We get 

2/d 

(99) 



2?r 


n 


U 










-01 


m 


2 r 


\1~ 





We determine €f for the interacting gas by inserting n from equ. (|97|l into 
equ. Ij99(l . The universal parameter is ^ = /u/ef • We find 

£ = ie 3/2 + ^e 5/2 log(e) - 0.025e 5 / 2 + . . . = 0.475 (e = 1), (100) 
2 16 

which agrees quite well with the result of fixed node quantum Monte Carlo cal- 
culations. The calculation has been extended to 0(e 7 / 2 ) by Arnold et al. |47j . 
Unfortunately, the next term is very large and it appears necessary to combine 
the expansion in 4 — e dimensions with a 2 + e expansion in order to extract 
useful results. The e expansion has also been applied to the calculation of 
the gap the critical temperature 0H] and the critical chemical potential 
imbalance |4^1l4^| . 
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5 QCD and its symmetries 
5.1 Introduction 

Before we discuss QCD at large baryon density we would like to provide a 
quick review of QCD and the symmetries of QCD. The elementary degrees 
of freedom are quark helds i\j a a j and gluons A 1 ^ . Here, a is color index that 
transforms in the fundamental representation for fermions and in the adjoint 
representation for gluons. Also, / labels the quark flavors u,d 7 s 7 c,b,t. In 
practice, we will focus on the three light flavors up, down and strange. The 
QCD lagrangian is 

C = ^ f {ip-m f )i> f --G^G^, (101) 
/ 

where the field strength tensor is defined by 

G% = d^Al - d v Al + gf abc AlAt, (102) 
and the covariant derivative acting on quark fields is 

m = ^[id^ + gA a ^j^. (103) 

QCD has a number of remarkable properties. Most remarkably, even though 
QCD accounts for the rich phenomenology of hadronic and nuclear physics, 
it is an essentially parameter free theory. To first approximation, the masses 
of the light quarks u, d, s are too small to be important, while the masses 
of the heavy quarks c, b, t are too heavy. If we set the masses of the light 
quarks to zero and take the masses of the heavy quarks to be infinite then 
the only parameter in the QCD lagrangian is the coupling constant, g. Once 
quantum corrections are taken into account g becomes a function of the scale 
at which it is measured. If the scale is large then the coupling is small, but in 
the infrared the coupling becomes large. This is the famous phenomenon of 
asymptotic freedom. Since the coupling depends on the scale the dimensionlcss 
parameter g is traded for a dimensionful scale parameter Aqcd- Since Aqcd 
is the only dimensionful quantity in QCD with massless fermions it is not 
really a parameter of QCD, but reflects our choice of units. In standard units, 
Aqcd ^ 200MeV~ lfm -1 . 

Another important feature of the QCD lagrangian are its symmetries. 
First of all, the lagrangian is invariant under local gauge transformations 
U(x) g SU(3) C 

iP{x) - U(x)i/>{x), A^x) - U{x)A ll U\x) + lU(x)d^uHx), (104) 

where A^ = A^(\ a /2). In the QCD ground state at zero temperature and 
density the local color symmetry is confined. This implies that all excitations 
are singlets under the gauge group. 
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The dynamics of QCD is completely independent of flavor. This implies 
that if the masses of the quarks are equal, m u = = m s , then the theory is 
invariant under arbitrary flavor rotations of the quark fields 

4>f "> VfgiPg, (105) 

where V £ SU(3). This is the well known flavor (isospin) symmetry of the 
strong interactions. If the quark masses are not just equal, but equal to zero, 
then the flavor symmetry is enlarged. This can be seen by defining left and 
right-handed fields 

^L,fl = i(l±75)V>. (106) 
In terms of L/R fields the fermionic lagrangian is 

C = $ L {%Q)1) L + i) R {iIp)i)R + 4> L M?p R + iPrMtPl, (107) 

where M = dia,g(m u ,md,m s ). We observe that if quarks are massless, 
m u = irid = m s = 0, then there is no coupling between left and right handed 
fields. As a consequence, the lagrangian is invariant under independent flavor 
transformations of the left and right handed fields. 

^lj -» L fg ip L , g , ip R j -» R fg ipR, g , (108) 

where (L,R) £ SU(3)l x SU(3)r. In the real world, of course, the masses of 
the up, down and strange quarks are not zero. Nevertheless, since m u , <C 
m s < Aqcd QCD has an approximate chiral symmetry. 

In the QCD ground state at zero temperature and density the flavor sym- 
metry is realized, but the chiral symmetry is spontaneously broken by a quark- 
anti-quark condensate {'<PliPr + '4>riPl)- As a result, the observed hadrons can 
be approximately assigned to representations of the SU(3)v flavor group, but 
not to representations of SU(3)l x SU(3)r. Nevertheless, chiral symmetry 
has important implications for the dynamics of QCD at low energy. Gold- 
stone's theorem implies that the breaking of SU(3)l x SU(3)r — > SU(3)v is 
associated with the appearance of an octet of (approximately) massless pseu- 
doscalar Goldstone bosons. Chiral symmetry places important restrictions on 
the interaction of the Goldstone bosons. These constraints are obtained most 
easily from the low energy effective chiral lagrangian. At leading order we have 

C = ^LTr [d^Sd^S^] + |BTr(Mi;t) + h.cj + ..., (109) 

where £ — exp(z</< a A a // 7r ) is the chiral field, f n is the pion decay constant 
and M is the mass matrix. Expanding S in powers of the pion, kaon and eta 
fields 4> a we can derive the leading order chiral perturbation theory results for 
Goldstone boson scattering and the coupling of Goldstone bosons to external 
fields. Higher order corrections originate from loops and higher order terms 
in the effective lagrangian. 
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Finally, we observe that the QCD lagrangian has two U(l) symmetries, 

U(1) B - i>L - e^ L , tpR — > e^ipR (110) 

U(1) A : ^L^e M Vi, ^«-»e- ia Vfl. (HI) 

The J7(l)s symmetry is exact even if the quarks are not massless. Superficially 
it appears that the U(1)a symmetry is explicitly broken by the quark masses 
and spontaneously broken by the quark condensate. However, there is no 
Goldstone boson associated with spontaneous U(1)a breaking. The reason is 
that at the quantum level the U(1)a symmetry is broken by an anomaly. The 
divergence of the U(1)a current is given by 

- S? G ^, (112) 

where G£„ = e A1 ^ a( gG^ /3 /2 is the dual field strength tensor. 



5.2 QCD at finite density 

In the real world the quark masses are not equal and the only exact global 
symmetries of QCD are the U(l) / flavor symmetries associated with the con- 
servation of the number of up, down, and strange quarks. If we take into 
account the weak interactions then flavor is no longer conserved and the only 
exact symmetries are the U(1)b of baryon number and the U(1)q of electric 
charge. 

In the following we study hadronic matter at non-zero baryon density. 
We will mostly focus on systems at non-zero baryon chemical potential but 
zero electron U(1)q chemical potential. We should note that in the context 
of neutron stars we are interested in situations when the electric charge, but 
not necessarily the electron chemical potential, is zero. Also, if the system 
is in equilibrium with respect to strong, but not to weak interactions, then 
non-zero flavor chemical potentials may come into play. 

The partition function of QCD at non-zero baryon chemical potential is 
given by 

z = y> P (-«), 

where i labels all quantum states of the system, Ei and Ni are the energy and 
baryon number of the state i. If the temperature and chemical potential are 
both zero then only the ground state contributes to the partition function. 
All other states give contributions that are exponentially small if the volume 
of the system is taken to infinity. In QCD there is a mass gap for states that 
carry baryon number. As a consequence there is an onset chemical potential 



[i c = min(Ei/Ni), 



(114) 
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such that the partition function is independent of p for p < p c . For p > p c the 
baryon density is non-zero. If the chemical potential is just above the onset 
chemical potential we can describe QCD, to first approximation, as a dilute 
gas of non- interacting nucleons. In this approximation p c — mpf. Of course, 
the interaction between nucleons cannot be neglected. Without it, we would 
not have stable nuclei. As a consequence, nuclear matter is self-bound and the 
energy per baryon in the ground state is given by 

~~j^r — m N — — 15 MeV. (115) 

The onset transition is a first order transition at which the baryon density 
jumps from zero to nuclear matter saturation density, po — 0.14fm~ 3 . The 
first order transition continues into the finite temperature plane and ends at 
a critical endpoint at T = T c ~ 10 MeV, see Fig. El 

Nuclear matter is a complicated many-body system and, unlike the situ- 
ation at zero density and finite temperature, there is little information from 
numerical simulations on the lattice. This is related to the so-called 'sign 
problem'. At non-zero chemical potential the euclidean fermion determinant 
is complex and standard Monte-Carlo techniques based on importance sam- 
pling fail. Recently, some progress has been made in simulating QCD for small 
p and T ~ T c [501 1511 152j , but the regime of small temperature remains inac- 
cessible. 

In neutron stars there is a non-zero electron chemical potential and matter 
is neutron rich. Pure neutron matter has positive pressure and is stable at 
arbitrarily low density. As we emphasized in Sect.0]dilute neutron matter has 
universal properties that can be explored using atomic systems. As the density 
increases three and four-body interactions as well as short range forces become 
more important and effective field theory methods are no longer applicable. 

If the density is much larger than nuclear matter saturation density, p 3> 
Po, we expect the problem to simplify. In this regime it is natural to use 
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a system of non-interacting quarks as a starting point |53j . The low energy 
degrees of freedom are quark excitations and holes in the vicinity of the Fermi 
surface. Since the Fermi momentum is large, asymptotic freedom implies that 
the interaction between quasi-particles is weak. We shall see that this does 
not imply that the phase diagram is simple, but it does imply that the phase 
structure can be studied in a systematic fashion. 

6 Effective field theory near the Fermi surface 
6.1 High density effective theory 

The QCD Lagrangian in the presence of a chemical potential is given by 

C = $ (tip + M7 o ~ M) i> - \G;„G a ^, (116) 

where = + igA^ is the covariant derivative, M is the mass matrix 
and fj, is the baryon chemical potential. If the baryon chemical potential is 
large, /x 3> Aqcd, then we expect the effective coupling to be small and 
weak coupling methods to be applicable. We shall see, however, that the weak 
coupling expansion is not a simple expansion in the number of loops. Effective 
field theory methods are useful in constructing a systematic weak coupling 
expansion. 

The main observation is that the relevant low energy degrees of freedom 
are particle and hole excitations in the vicinity of the Fermi surface. We shall 
describe these excitations in terms of the field ip v (x), where v is the Fermi 
velocity. At tree level, the quark field ip can be decomposed as ip = ip Vi +-\-i() Vi - 
where ip v> ± = P v ,±ip with P v ,± — ^(l±ot-v)ip. Note that P v ,± is a projector on 
states with positive/negative energy. To leading order in l//twe can eliminate 
the field ip- using its equation of motion. The lagrangian for the ip+ field is 
given by [SI EH EE] 

C = ^^ v .D-^-^^ + ..^-\g^G^ + .... (117) 

with = (1,1?). Note that v labels patches on the Fermi surface, and that 
the number of these patches grows as /x 2 . The leading order v ■ D interaction 
does not connect quarks with different v, but soft gluons can be exchanged 
between quarks in different patches. In addition to that, there are four, six, 
etc. fermion operators that contain fermion fields with different velocity labels. 
These operators are constrained by the condition that the sum of the velocities 
has to be zero. 

In the case of four-fermion operators there are two kinds of interactions 
that satisfy this constraint, see Fig. The first possibility is that both the 
incoming and outgoing fermion momenta are back-to-back. This corresponds 
to the BCS interaction 
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Fig. 10. Kinematics of four-fermion operators in the effective theory. 

c = ^ J2 vr'Rr\v.v>)^ v r^ v )(i,t,r'^_ v ,), (us) 

v',r,r' 

where v ■ v' = cos 9 is the scattering angle, R[ r (x) is a set of orthogonal 
polynomials, and r, r' determine the color, flavor and spin structure. The 
second possibility is that the final momenta are equal to the initial momenta 
up to a rotation around the axis defined by the sum of the incoming momenta. 
The relevant four-fermion operator is 

c = - 2 E Ff r 'R{ r '(v ■ t/)(y„A/v) (4r'4). (H9) 

v',r,r' 

In a system with short range interactions only the quantities -Fz(O) are known 
as Fermi liquid parameters. 

6.2 Hard Loops 

The effective field theory expansion is complicated by the fact that the number 
of patches N v ~ n 2 /A 2 grows with the chemical potential. This implies that 
some higher order contributions that are suppressed by l//i 2 can be enhanced 
by powers of N v . The natural solution to this problem is to sum the leading 
order diagrams in the large N v limit |57j . For gluon n-point functions this 
corresponds to the well known hard dense loop approximation 58 59 60 . 

The simplest example is the gluon two point function. At leading order in 
g and we have 

II? v (p) = 2g 2 N f 6 -^-J2 v » v » I 7?T47I FW7 r~V ^ 

2 ^ J (27r) 4 (fco - tfc)(«o +Po - h+p) 

where I). = v ■ k. We note that taking the momentum of the external gluon 
to zero automatically selects forward scattering. We also observe that the 
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Fig. 11. Hard dense loop contribution to gluon n-point functions. 

gluon can interact with fermions of any Fermi velocity so that the polarization 
function involves a sum over all patches. After performing the ko integration 
we get 



rafc 



np{p) = 2g*N f 



dh 



In 



dn k 



(2tt) 2 J 2?r po 



l p dlk 



(121) 



where n k is the Fermi distribution function. We note that the l k integration is 
automatically restricted to small momenta. The integral over the transverse 
momenta l±, on the other hand, diverges quadratically with the cutoff A±. 
We observe, however, that the sum over patches and the integral over l± can 
be combined into an integral over the entire Fermi surface 



-E 

2tt ^ 



d 2 l ± 



JJ_ 

2tt 2 



df2 
An ' 



(122) 



This means that the transverse momentum integral is extended all the way 
up to fi. Because the energy of the fermions is small but the loop momentum 
is large the integral is referred to as a hard dense loop. We find 



nf v (p) = 2m V / ^v,v v 



dn 

An ' 



1 - 



Po 



Po 



I, 



(123) 



where we have defined the effective gluon mass m 2 — Np g 2 fj, 2 / (Att 2 ) . This re- 
sult has the same structure as the non-relativistic expression given in equ. (|33|l . 
but the tadpole contribution is missing. As a consequence, equ. (|123fl is not 
transverse. In the relativistic theory the tadpole contribution originates from 
the D 2 _/(2fi) in the effective lagrangian. The tadpole is proportional to the 
total density and corresponds to a counterterm |54j 



„ 1 9 f dQ . . 



(124) 



Putting everything together we find 



n*l(p)=2m 2 S ab J |J{<WU 



Po - lp 



(125) 
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The gluonic three-point function shown in Fig. Illb can be computed in the 
same fashion. We find 

(126) 

where p,q,r are the incoming gluon momenta (p + q + r = 0) . We note that in 
the case of the three point function, as well as in all higher n-point functions, 
there are no tadpole or counterterm contributions. There is a simple generat- 
ing functional for these loop integrals which is known as the hard dense loop 
(HDL) effective action |rjl| 

£hdl = G U (v. V D )2 G rf- ( 12? ) 

This is a gauge invariant, but non-local, effective lagrangian. 



6.3 Non-Fermi liquid effective field theory 

In this Section we shall study the effective field theory in the regime oj < m 
where u> is the excitation energy and m is the effective gluon mass |fi2j . In 
the previous section we argued that hard dense loops have to be resummed 
in order to obtain a consistent low energy expansion. The effective lagrangian 
is given by 

C = V4 (iv-D- ^ - l -G%G% + C HDL + £ 4/ + . . . . (128) 

Since electric fields are screened the interaction at low energies is dominated 
by the exchange of magnetic gluons. The transverse gauge boson propagator 
is 

A,(fc) = - „ ^o-hki) 

where we have assumed that \ko\ < \k\. We observe that the propagator 
becomes large in the regime |fcrj| ~ \k\ 3 /m 2 . If the energy is small, |fcrj| "C m, 
then the typical energy is much smaller than the typical momentum, 

|fe|^(m 2 |fc |) 1 / 3 »|A :o |. (130) 

This implies that the gluon is very far off its energy shell and not a propagating 
state. We can compute loop diagrams containing quarks and transverse gluons 
by picking up the pole in the quark propagator, and then integrate over the 
cut in the gluon propagator using the kinematics dictated by equ. (|130|l . In 
order for a quark to absorb the large momentum carried by a gluon and stay 
close to the Fermi surface the gluon momentum has to be transverse to the 
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Fig. 12. One-loop contributions to the quark self energy and the quark-gluon vertex. 
The black blob in the third diagram denotes the HDL gluon three point function. 
In the magnetic regime the graphs scale as ojlog(oj), ui 1 ^ 3 and w 2 ^ 3 , respectively. 

momentum of the quark. This means that the term fcj_/(2/i) in the quark 
propagator is relevant and has to be kept at leading order. Equation l|13()(l 
shows that kj_/(2(j,) ^> kg as ko — > 0. This means that the pole of the quark 
propagator is governed by the condition kn ~ fc^/(2/i). We find 

k^ ff 2 V /3 fco /3 , fc„~ 5 4/ V /3 fco /3 - (131) 

In the low energy regime propagators and vertices can be simplified even 
further. The quark and gluon propagators are 

Safi(p) = ^ , (132) 

Po-P\\ - 2^ +iesgn(p a ) 

Aj(fc) = -5 — — e-, (133) 

fci-ifm 2 It 

and the quark gluon vertex is gvi(\ a /2). Higher order corrections can be found 
by expanding the quark and gluon propagators as well as the HDL vertices in 
powers of the small parameter e = (ko/m). 

The regime in which all momenta, including external ones, satisfy the 
scaling relation is completely perturbative, i.e. graphs with extra loops 

are always suppressed by extra powers of e 1 / 3 . One way to see this is to 
rescale the fields in the effective lagrangian so that the kinetic terms are scale 
invariant under the transformation (xq,xu,x±) — ► (e _1 £o> e 2 / 3 £C| | , £ _1 / 3 x_l). 
The scaling behavior of the fields is ^ — ^ e 5 ^ 6 ?A and At — > e 5 ^ 6 Ai. We find 
that the scaling dimension of all interaction terms is positive. The quark 
gluon vertex scales as e 1 / 6 , the HDL three gluon vertex scales as e 1 / 2 , and the 
four gluon vertex scales as e. Since higher order diagrams involve at least one 
pair of quark gluon vertices the expansion involves positive powers of e 1 / 3 . 

As a simple example we consider the fermion self energy in the limit po — > 
0. The one-loop diagram is 
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. 2 f dk f dkj_ kj_ 

dk \\ O(p + k ) 



2vr fe|| +p|| - (k x + p 1 _) 2 /{2^i) +»e' 



(134) 



with Ci? = (N% — 1)/(2N C ) and 77 = (7r/2)m 2 . This expression shows a number 
of interesting features. First we observe that the longitudinal and transverse 
momentum integrations factorize. The longitudinal momentum integral can 
be performed by picking up the pole in the quark propagator. The result 
is independent of the external momenta and only depends on the external 
energy. The transverse momentum integral is logarithmically divergent. We 

find [HHI EH ESI ESI EZ| 

27 (p) = #aPologfAV (135) 



where A is a cutoff for the k± integral. We showed that the logarithmic diver- 
gence can be absorbed in the parameters of the effective theory [37| . In order 
to fix the corresponding counterterm we have include electric gluon exchanges 
For fco <C m the electric gluon propagator is given by 

Do,(k) = -pr^?. (136) 

Higher order corrections can be obtained by expanding the full HDL expres- 
sion in powers of ko/m. The electric contribution is dominated by large mo- 
menta and does not contribute to fractional powers or logarithms of ko- We 
get 

dkn f dk 2 , 1 



S(p) =ig C F ' ' 



2tt J (2tt) 2 k 2 _ + 2m 2 
dk\\ O(p + k ) 



2tt fey - (kx +pi_) 2 /(2fi) +ie 



(137) 



This contribution scales as po\og(A/m). The logarithm of the cutoff cancels 
the logarithm in equ. ^ 135(1 . We get [S7| 



. ( 4\/2m\ .it. 

Polog [^\M +P0 + t 2 M 



(138) 



Finally, there are contributions from the hard regime in which both the en- 
ergy and the momentum of the gluon are large, ko ~ \k\ > m. This regime 
corresponds to the HDL term in the fermion self energy EE] ■ The HDL 
term gives an 0(g 2 ) correction to the low energy parameters Vf and Sfi. 

The logarithmic term in the fermion self energy leads to a breakdown of 
Fermi liquid theory. The quasi-particle velocity 
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Fig. 13. Gap equation for the superfluid gap in a theory with long range interac- 
tions. 

u(po) = TT^M (139) 

and the wave function renormalization go to zero logarithmically as the quasi- 
particle energy goes to zero. One physical consequence of this behavior is an 
anomalous Tlog(T) term in the specific heat j69H67j . The effective theory can 
also be used to study perturbative corrections in other quantities. We find, 
in particular, a QCD version of Migdal's theorem. Migdal showed that vertex 
corrections to the electron-phonon coupling are suppressed by the ratio of the 
electron mass to the mass of the positive ions 0]. In the Landau damping 
regime of QCD loop corrections to the quark-gluon vertex are suppressed by 
powers of e 1 / 3 . 



6.4 Color superconductivity 

In Sect. H3.1(l we showed that the particle-particle scattering amplitude in the 
BCS channel q(p) + q{—p) — * q{p') + q(—p') is special. The total momentum 
of the pair vanishes and as a consequence loop corrections to the scattering 
amplitude are enhanced. This implies that all ladder diagrams have to be 
summed. Crossed ladders, vertex corrections, etc. involve momenta in the 
regime k± k\\ ^ ko and are perturbative. 

If the interaction in the particle-particle channel is attractive then the BCS 
singularity leads to the formation of a pair condensate and to a gap in the 
fcrmion spectrum. The gap can be computed by solving a Dyson-Schwinger 
equation for the anomalous (particle-particle) self energy. In QCD the inter- 
action is attractive in the color anti-triplet channel. The structure of the gap 
is simplest in the case of two flavors. In that case, there is a unique color 
anti-symmetric spin zero gap term of the form 

(CC75^)=^e,, (140) 

Here, a, b labels color and i, j flavor. The gap equation is given by 

a f x • 9 ^> f dk f dk 2 : kj_ 
A(p ) = -ig 2 C A 1 ' 



2vr J (2tt) 2 k 3 ± + r)(k - pa) 

(Ml) 

2tt k 2 + k 2 +A(k ) 2 ' [ ' 
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where Ca = 2/3 is a color factor. Like the normal self energy, the anomalous 
self energy A(p) is dominantly a function of energy. We carry out the integrals 
over k± and ku and analytically continue to imaginary energy p 4 = ip a . The 
euclidean gap equation is [70] [7J [73 EH 

A{pA) = Jl_ f dk4 l 0g ( A BCS \ *M ( } 

The scale Abcs is sensitive to electric gluon exchange. In the anomalous 
self energy the logarithmic divergence does not cancel between magnetic 
and electric gluon exchanges. The reason is that the magnetic contribu- 
tion is proportional to SijViVj in the normal self energy and 5ijVi(—Vj) in 
the anomalous case. The logarithmic dependence on the cutoff is absorbed 
by the BCS four-fermion operator. A simple matching calculation gives 
Abcs = 2567r 4 (2/A r /) 5/2 g~ f > :68:. The solution to the gap equation was 
found by Son [70] 

A M = A°Hak s ) (W3) 

where x = log(2ABcs /(p4 + w p ) and w 2 = p\ + Aq. The gap on the Fermi 
surface is 

4\ ( 3tt 2 \ 



A q ~ 2 Abcs exp ^ — J exp y-~^J ■ ( 144 ) 

This result is correct up to 0{g) corrections to the pre-exponent. In order 
to achieve this accuracy the g 2 po log(po) term in the normal self energy, 
equ. I|135|l , has to be included in the gap equation [7^ [TSJ [HS] . 

The order parameter is slightly more complicated in QCD with Nf = 3 
massless flavors. The energetically preferred phase is the color-flavor-locked 
(CFL) phase described by [2^] 

«C<yfe$) = * (Wj ~ S ! S T) ■ (145) 

In the CFL phase there are eight fermions with gap Aqfl and one fermion 
with gap 2A CFL . The CFL gap is given by A C fl = 2" 1/3 Z\ [77j. The CFL 
phase has a number of remarkable properties [7S1 [7S] . Most notably, chiral 
symmetry is broken in the CFL phase and the low energy spectrum contains 
a flavor octet of pseudoscalar bosons. We shall study the dynamics of these 
Goldstone modes in Sect.[7J 



6.5 Mass terms 



Mass terms modify the parameters in the effective lagrangian. These parame- 
ters include the Fermi velocity, the effective chemical potential, the screening 
mass, the BCS terms and the Landau parameters. At tree level the correction 
to the Fermi velocity and the chemical potential are given by 
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Fig. 14. Mass terms in the high density effective theory. The first diagram shows a 
O(MM^) term that arises from integrating out the tp- field in the QCD lagrangian. 
The second diagram shows a 0(M 2 ) four-fermion operator which arises from inte- 
grating out ip- and hard gluon exchanges. 



«*■ = !-£*-. ^ = ~- (146) 

%P F 2pf 

The shift in the Fermi velocity also affects the coupling gvF of a magnetic 
gluon to quarks. It is important to note that at leading order in g this the 
only mass correction to the coupling. This is not entirely obvious, as one can 
imagine a process in which the quark emits a gluon, makes a transition to 
a virtual high energy state, and then couples back to a low energy state by 
a mass insertion. This process would give an 0(m//i) correction to g, but it 
vanishes in the forward direction |79| . 

Quark masses modify quark-quark scattering amplitudes and the corre- 
sponding Landau and BCS type four-fermion operators. Consider quark-quark 
scattering in the forward direction, v + v' — > v + v' . At tree level in QCD this 
process receives contribution from the direct and exchange graph. In the ef- 
fective theory the direct term is reproduced by the collinear interaction while 
the exchange terms has to be matched against a four-fermion operator. The 
spin-color-flavor symmetric part of the exchange amplitude is given by 



M(v, v>; v, v>) = \i-'^-±-\ (147) 



Cf 9 2 f m 2 
AN c N f p 2 F \ p%l 

where Cf — (N^ — V)/(2N C ) and x = v ■ v' is the scattering angle. We observe 
that the amplitude is independent of x in the limit m — > 0. Mass corrections 
are singular as x — > 1. The means that the Landau coefficients Fi contain 
logarithms of the cutoff. We note that there is one linear combination of 
Landau coefficients, Fo — i*i/3, which is cutoff independent. 

Equations 1 1461 147(1 are valid for Nf > 1 degenerate flavors. Spin and 
color anti-symmetric BCS amplitudes require at least two different flavors. 
Consider BCS scattering v + (—v) — > v' + (—v 1 ) in the helicity flip channel 
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L + L — » R + R. The color-anti-triplet amplitude is given by 

M(v,-v;v'-v') = (148) 

^ Pp I J F 

where mi and ni2 are the masses of the two quarks and Ca = [N c + 1 )/ (2N C ) . 
We observe that the scattering amplitude is independent of the scattering 
angle. This means that at leading order in g and m only the s-wave potential 
Vo is non-zero. 

In order to match Green functions in the high density effective theory to 
an effective chiral theory of the CFL phase we need to generalize our results 
to a complex mass matrix of the form C = —^lMi/jji — ipj^M^ip^, see Fig. 1141 
The dfi term is 

C = ~ {^l+ mm ^l+ + ^ R+ M^M^ R+ ) . (149) 
and the four-fermion operator in the BCS channel is 

c = Kr c <r) [W ac W bd {M) lk {M) 3l ] 

+ (L<r+R,M <-> M f ) . (150) 

7 Chiral theory of the CFL phase 

7.1 Introduction 

For energies smaller than the gap the only relevant degrees of freedom are the 
Goldstone modes associated with spontaneously broken global symmetries. 
The quantum numbers of the Goldstone modes depend on the symmetries 
of the order parameter. In the following we shall concentrate on the CFL 
phase. Goldstone modes determine the specific heat, transport properties, 
and the response to external fields for temperatures less than T c . As we shall 
see, Goldstone modes also determine the phase structure as a function of the 
quark masses. 

7.2 Chiral effective field theory 

In the CFL phase the pattern of chiral symmetry breaking is identical to the 
one at T = fj, = 0. This implies that the effective lagrangian has the same 
structure as chiral perturbation theory. The main difference is that Lorentz- 
invariance is broken and only rotational invariance is a good symmetry. The 
effective lagrangian for the Goldstone modes is given by [201 
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f 2 

C ef/ = ^Tr [WoSWoS^ - vldiSdiS^] + [BTr(M^) + h.c] 

+ [AiTr(MX ,t )Ti(MZ ,t ) + A 2 Tr(MS^MS^) 

+A 3 Tr(Mi: t )Tr(M t S) + h.c] + .... (151) 

Here E = exp(i^> a A a // w ) is the chiral field, /„• is the pion decay constant and 
M is a complex mass matrix. The chiral field and the mass matrix transform as 
S — > LJCi? 1 ' and M — » LMB) under chiral transformations (L, i?) S SU(3) L x 
5£/(3)fl. We have suppressed the singlet fields associated with the breaking 
of the exact U(l)v and approximate U(1)a symmetries. 

At low density the coefficients f w , B, A,, . .. are non-perturbative quan- 
tities that have to extracted from experiment or measured on the lattice. 
At large density, on the other hand, the chiral coefficients can be calculated 
in perturbative QCD. The pion decay constant and the pion velocity can be 
determined by gauging the ST/ (3) l x SU(3)r symmetry. The covariant deriva- 
tive D^S = d^S + iW^S — i£W^ generates mass terms for the gauge field 
W L ' R 

c = T G W - w & - T W ~ w ^) ■ (152) 

The electric and magnetic screening masses in the CFL phase can be deter- 
mined as in Sect. 13.41 The main difference is that in the CFL phase there are 
nine different fermion modes, and that not all of these modes have the same 
gap. There is also mixing between flavor and color gauge fields. It is easiest to 
compute the screening for the color gauge fields. The electric screening mass 
is 



n m = -2i 



d 4 p I 7 P 2 + e 2 p 
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3(p 2 -e 2 p -A 2 )(p 2 -e 2 p -4A 2 ) 



(153) 



The first terms comes from particle-hole diagrams with two octet quasi- 
particles while the second term comes from diagrams with one octet and 
one singlet quasi-particle. There is no coupling of an octet field to two singlet 
particles. The third and fourth term are the corresponding contributions from 
particle-particle and hole-hole pairs. In the CFL phase A\ — 2A% = 2A. The 
four integrals in (|153fl give 

'6 3 3 9 J Utt 2 J 18 V27T 2 ' 



The magnetic mass can be computed in the same fashion. As in Sect. 13.41 we 
have to add the contribution from the tadpole and the structure of the gauge 
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Fig. 15. Contribution of the 0{M 2 ) BCS four-fermion operator to the condensation 
energy in the CFL phase . 



field mass term is m 2 D {A^ — A 2 /3). Mixing between flavor and color gauge 



fields was studied in jSOJ |^| . The result is that there is no screening for the 
vector field Wl + Wr and that the screening mass for the axial field Wl — Wr 
is equal to the mass of the color gauge field. We conclude that 
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(155) 



Mass terms are determined by the operators studied in Sect. 16.51 We observe 
that both equ. Ijl49(l and (|150|l are quadratic in M. This implies that B = in 
perturbative QCD. B receives non-perturbative contributions from instantons, 
but these effects are small if the density is large, see |83) . 

We also note that X L = MM^ /(2p F ) and X R = M^M/(2p F ) in equ. I|149|) 
act as effective chemical potentials for left and right-handed fermions, re- 
spectively. Formally, the effective lagrangian has an SU(3)l x SU(3)r gauge 
symmetry under which r transform as the temporal components of non- 
abelian gauge fields. We can implement this approximate gauge symmetry in 
the CFL chiral theory by promoting time derivatives to covariant derivatives 

ED, 



/MMt\ 



(156) 



The BCS four-fermion operator in equ. (|150|l contributes to to the condensa- 
tion energy in the CFL phase, see Fig. El The diagram is proportional to the 
square of the superfluid density 



5*81 
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(2tt) 4 p 2 



A 2 ( Pa ) 



JJ_ 

2tt 2 



(157) 



We note that the superfluid density is sensitive to energies p > A and the 
energy dependence of the gap has to be kept. The color- favor factor is 
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= -l{(Tr[M]) 2 -Tr 



(158) 



We also note that the four-fermion operator is proportional to g 2 and the 
explicit dependence of the diagram on g cancels. We find E5] 



AS 



3A 2 
'4^2" 



(Tr[M])' 



Tr 



M 



r2 



+ (m <-> M f ) . (159) 



This term can be matched against the Ai terms in the effective lagrangian. 
The result is |521I75| 



-An 



3A 2 



^3 = 0. 



(160) 



We can now summarize the structure of the chiral expansion in the CFL 
phase. The effective lagrangian has the form 



C 



(161) 



Loop graphs in the effective theory are suppressed by powers of 9/(47r/ ff ). 
Since the pion decay constant scales as f„ ~ pp Goldstone boson loops are 
suppressed compared to higher order contact terms. We also note that the 
quark mass expansion is controlled by m 2 / (ppA). This is means that the 
chiral expansion breaks down if m 2 ~ pfA. This is the same scale at which 
BCS calculations find a transition from the CFL phase to a less symmetric 
state. 



7.3 Kaon condensation 

Using the chiral effective lagrangian we can determine the dependence of the 
order parameter on the quark masses. We will focus on the physically relevant 
case m s > m u = md- Because the main expansion parameter is m 2 /(pFA) 
increasing the quark mass is roughly equivalent to lowering the density. The 
effective potential for the order parameter is 



V eff = ^Tr [2X L SX R S^ - X 2 - X 2 R ]-A 1 [(Tr(M&)) 2 - Tr ((MS^f 

(162) 

The first term contains the effective chemical potential /j, s — m 2 /(2pp) and 
favors states with a deficit of strange quarks (with strangeness S = —1). The 
second term favors the neutral ground state S = 1. The lightest excitation 
with positive strangeness is the K° meson. We therefore consider the ansatz 
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Fig. 16. Phase structure of CFL matter as a function of the effective chemical 
potential /i s = m 2 s /{2pF) and the lepton chemical potential fj,Q, from Kaplan & 
Reddy (2001). A typical value of fj, s in a neutron star is 10 MeV. 



S = exp(iaA4) which allows the order parameter to rotate in the K° direction. 
The vacuum energy is 



V(a) =-%[- 



1 f mi 



2pi 



sin(a) 2 + {m° K ) 2 (cos(a) - 1) , (163) 



where (mP K ) 2 = (AAi/ f 2 )m{m + m s ). Minimizing the vacuum energy we ob- 
tain 

( 1 n 8 < m° K 



cos(a) 

The hypercharge density is 



/i, > TO 



(164) 



K 



«) 4 

Ml 



(165) 



This result has the same structure as the charge density of a weakly interacting 
Bose condensate. Using the perturbative result for A\ we can get an estimate 
of the critical strange quark mass. We find 



m s {crit) = 3.03 • m^A 2 ^, 



(166) 



from which we obtain m s (crit) ~ 70 MeV for A ~ 50 MeV. This result sug- 
gests that strange quark matter at densities that can be achieved in neutron 
stars is kaon condensed. We also note that the difference in condensation en- 
ergy between the CFL phase and the kaon condensed state is not necessarily 
small. For fj, s — > A we have sin(a) — > 1 and V(a) — » f 2 A 2 /2. Since f 2 is 
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Fig. 17. This figure shows the fermion spectrum in the CFL phase. For m 3 — there 
are eight fermions with gap A and one fermion with gap 2 A (not shown). Without 
kaon condensation gapless fermion modes appear at fi s = A (dashed lines). With 
kaon condensation gapless modes appear at fi s — 4A/3. 

of order /i 2 /(27r 2 ) this is comparable to the condensation energy in the CFL 
phase. 

The strange quark mass breaks the SU(3) flavor symmetry to SU(2)i x 
U(1)y- In the kaon condensed phase this symmetry is spontaneously broken 
to U(1)q. If isospin is an exact symmetry there are two exactly massless 
Goldstone modes [HI], the K° and the K + . Isospin breaking leads to a small 
mass for the K + . The phase structure as a function of the strange quark mass 
and non-zero lepton chemical potentials was studied by Kaplan and Reddy 
|85] , see Fig. ^| We observe that if the lepton chemical potential is non-zero 
charged kaon and pion condensates are also possible. 

7.4 Fermions in the CFL phase 

So far we have only studied Goldstone modes in the CFL phase. However, as 
the strange quark mass is increased it is possible that some of the fermion 
modes become light or even gapless [SS]. In order to study this question we 
have to include fermions in the effective field theory. The effective lagrangian 
for fermions in the CFL phase is |£3 |S5] 



C = Tr (NUv^D^N) - DTt (N^v^s N}) - FTt (7VV 75 [A^N]) 



Nl.r are left and right handed baryon fields in the adjoint representation of 
flavor SU(3). The baryon fields originate from quark-hadron complementarity 




(167) 
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|78j . We can think of N as describing a quark which is surrounded by a diquark 
cloud, N L ~ <Zl(<Zl<Zl)- The covariant derivative of the nucleon field is given 
by D^N = d^N + i[V M , N}. The vector and axial-vector currents are 

V^-^U^+C 1 ^}, 4, = (V,^) f, (168) 

where £ is defined by £ 2 = S. It follows that £ transforms as £ — > L^U(xy — 
U(x)t;R) with C/(x) S SC/^y- For pure SU(3) flavor transformations L = 
i? = V we have E/(x) = V. i* 1 and D are low energy constants that determine 
the baryon axial coupling. In perturbative QCD we find D = F = 1/2. 

The effective theory given in equ. (|167fl can be derived from QCD in the 
weak coupling limit. However, the structure of the theory is completely de- 
termined by chiral symmetry, even if the coupling is not weak. In particular, 
there are no free parameters in the baryon coupling to the vector current. 
Mass terms are also strongly constrained by chiral symmetry. The effective 
chemical potentials (Xl,Xr) appear as left and right-handed gauge potentials 
in the covariant derivative of the nucleon field. We have 

D N = d N + i[r ,N], (169) 
Fa = - l - {£ (do + iXn) ^ + & (do + iX L ) £} , 

where X L = MM^ l(2p F ) and X R = MHl/(2p F ) as before. (X L , X R ) covari- 
ant derivatives also appears in the axial vector current given in equ. 1168(1 . 

We can now study how the fermion spectrum depends on the quark mass. 
In the CFL state we have £ = 1. For fj, s — the baryon octet has an energy 
gap A and the singlet has gap 2A. The correction to this result comes from 
the term 

Tr (JVt[A„ N])) = ^ ((5-)Hz-) + (H°)t(E°) - (p)t(p) - (n)t(n)) , 

(170) 

where jl s — /i s diag(0, 0, 1). We observe that the excitation energy of the proton 
and neutron is lowered, ui p ^ n = A — fj, s , while the energy of the cascade states 
S~ , E!° particles is raised, u>~ = A + fi s . All other excitation energies are 
independent of fi s . As a consequence we find gapless (p, n) excitations at 
fi s = A. 

This result is also easily derived in microscopic models [SS]. The EFT 
perspective is nevertheless useful. In microscopic models the shift of the non- 
strange modes arises from a color chemical potential which is needed in or- 
der to neutralize the system. The effective theory is formulated directly in 
terms of gauge invariant variables and no color chemical potentials are needed. 
The shift in the non-strange modes is due to the fact that the gauge invari- 
ant fermion fields transform according to the adjoint representation of flavor 
517(3). 
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Fig. 18. Left panel: Energy density as a function of the current jk for several 
different values of ^t s = m 2 s /(2pF) close to the phase transition. Right panel: Ground 
state energy density as a function of fi s - We show the CFL phase, the kaon condensed 
CFL (KCFL) phase, and the supercurrent state (curKCFL). 

The situation is more complicated when kaon condensation is taken into 
account. In the kaon condensed phase there is mixing in the (p, S + , E~ , S~) 
and (n, E°, S°, A 8 , A ) sector. For <§C /i s <C A the spectrum is given by 



Numerical results for the eigenvalues are shown in Fig. El We observe that 
mixing within the charged and neutral baryon sectors leads to level repulsion. 
There are two modes that become light in the CFL window \i s < 2A. One 
mode is a linear combination of the proton and E + and the other mode is a 
linear combination of the neutral baryons (n, E°, S°, A 8 , A ). 

7.5 Meson supercurrent state 

Recently, several groups have shown that gapless fermion modes lead to insta- 
bilities in the current-current correlation function jHOllS]- Motivated by these 
results we have examined the stability of the kaon condensed phase against the 
formation of a non-zero current jHSIHSj- Consider a spatially varying U(l)y 
rotation of the maximal kaon condensate 




(171) 




1 

1/V2 i e ^K{x)j 

Q{ e -i<t>K{x) 1^/2 l/y/2 




(172) 



This state is characterized by non-zero currents 




(173) 
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In the following we compute the vacuum energy as a function of the kaon 
current jx = "^4>k- The meson part of the effective lagrangian gives a positive 
contribution 

£ = \vlfhl- (174) 

A negative contribution can arise from gapless fermions. In order to determine 
this contribution we have to calculate the fermion spectrum in the presence 
of a non-zero current. The relevant part of the effective lagrangian is 

C = Tr (NUv^D^N) + Tr (N^ 5 (p A + v ■ A) N) 

+ y {Tr (NN) - Tr (JV) Tr (N) + h.c.} , (175) 

where we have used D = F = 1/2. The covariant derivative is DqN = doN + 
i\p v , N] and DiN = d t N + iv ■ [V, N] with V, A given in equ. I JT75JI and 

^ = 2{^ t±et ^T { }- (176) 

The vector potential pv and the vector current V are diagonal in flavor space 
while the axial potential pa and the axial current A lead to mixing. The 
fermion spectrum is quite complicated. The dispersion relation of the lowest 
mode is approximately given by 

un^A+^^—\^-\v 3K , (177) 

where I = v -p — pp and we have expanded u>i near its minimum Iq = (p s +v ■ 
Jk)/4:. Equation l|177(l shows that there is a gapless mode if ^ s > 4A/3—jk/3. 
The contribution of the gapless mode to the vacuum energy is 



where dfl is an integral over the Fermi surface. The integral in equ. (1 178(1 re- 
ceives contributions from one of the pole caps on the Fermi surface. The result 
has exactly the same structure as the energy functional of a non-relativistic 
two-component Fermi liquid with non-zero polarization, see |94| . Introducing 
dimcnsionless variables 

x = —r, h= . (179) 

aA aA 

we can write E = cMfhix) with 

f h [x) =x 2 -- \{h + xf' 2 0{h +x)-{h- x) 5/2 0(h - x)] . (180) 



We have defined the constants 
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15 4 4 v % ' 



u 2 A 2 2 
^=^> a = ( 181 ) 



where = (21 — 81og(2))/36 is the numerical coefficient that appears in the 
weak coupling result for /„•. According to the analysis in |94| the function 
fh(x) develops a non-trivial minimum if hi < h < hi with hi ~ —0.067 and 
hi ~ 0.502. In perturbation theory we find a = 0.43 and the kaon condensed 
ground state becomes unstable for (A — 3/x s /4) < 0.007 A 

The energy density as a function of the current and the groundstate energy 
density as a function of fi s are shown in Fig. In these plots we have 
included the contribution of a baryon current jb, as suggested in |93| . In 
this case we have to minimize the energy with respect to two currents. The 
solution is of the form jb ~ 3k- The figure shows the dependence on jk for 
the optimum value oi jb- We have not properly implemented electric charge 
neutrality. Since the gapless mode is charged, enforcing electric neutrality will 
significantly suppress the magnitude of the current. We have also not included 
the possibility that the neutral mode becomes gapless. This will happen at 
somewhat larger values of /i s . 

We note that the ground state has no net current. This is clear from 
the fact that the ground state satisfies (5£/<5(V</>if ) = 0. As a consequence 
the meson current is canceled by an equal but opposite contribution from 
gapless fermions. We also expect that the ground state has no chromomagnetic 
instabilities. The kaon current is equivalent to an external gauge field. By 
minimizing the thermodynamic potential with respect to j we ensure that the 
second derivative, and therefore the screening mass, is positive. 



8 Conclusion: The many uses of effective field theory 

Strongly correlated quantum many body systems play a role in many different 
branches of physics, atomic physics, condensed matter physics, nuclear and 
particle physics. One of the main themes of these lectures is the idea that 
effective field theories provide a unified description of systems that involve 
vastly different scales. For example, nuclear physicists studying neutron mat- 
ter have learned a great deal from studying cold atomic gases (and vice versa). 
Similarly progress in understanding non-Fermi liquid behavior in strongly 
correlated electronic systems has been helpful in understanding dense quark 
matter in QCD. It is our hope that these lecture notes will play a small part 
in fostering exchange of ideas between different communities in the future. 
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